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The ultraperipheral heavy ion collisions (UPC's) are an important alternative to study 
the QCD dynamics until the next generation of e + e — /ep/eA colliders become reality. 
Due to the coherent action of all the protons in the nucleus, the electromagnetic field is 
very strong and the resulting flux of equivalent photons is large, which allows to study 
two-photon as well as photonuclcar interactions at high energies. In this paper we present 
a brief review of the vector meson production in UPC's at high energies using the QCD 
color dipole approach to describe their photonuclear production and the perturbative 
QCD Pomeron (BFKL dynamics) to describe the double meson production in photon- 
photon process. Predictions for rates and integrated cross sections are presented for 
energies of RHIC and LHC. 
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1. Introduction 

In ultraperipheral relativistic heavy-ion collisions (UPC's) the ions do not interact 
directly with each other and move essentially undisturbed along the beam direction. 
The only possible interaction is due to the long range electromagnetic interaction 
and diffractive processes (For a review see, e. g. Refs. ^El). Due to the coherent 
action of all the protons in the nucleus, the electromagnetic field is very strong 
and the resulting flux of equivalent photons is large. A photon stemming from the 
electromagnetic field of one of the two colliding nuclei can penetrate into the other 
nucleus and interact with one or more of its hadrons, giving rise to photon-nucleus 
collisions to an energy region hitherto unexplored experimentally. For example, 
the interaction of quasi-real photons with protons has been studied extensively at 
the electron-proton collider at HERA, with W lp < 200 GeV. Due to the larger 
number of photons coming from one of the colliding nuclei in heavy ion collisions, 
a similar and more detailed study will be possible in these collisions, with Wjn 



2 V. P. Gongalves and M. V. T. Machado 



reaching 950 GeV for the Large Hadron Collider (LHC) operating in its heavy 

ion mode. Similarly, one can inv estiga tes QCD effects in the context of coherent 

two-photon scattering at UPC's The important advantage of investigating 

two-photon interactions in nuclear collisions is the flux of virtual photons from the 

electromagnetic field of the nuclei scaling as Z 2 , giving a two-photon luminosity 

scaling as Z 4 . The maximum 77 collision energy W 77 is 2j/Ra, about 6 GeV at 

RHIC and 200 GeV at LHC 03 l n particular, the LHC will have a significant energy 

and luminosity reaching beyond LEP2, and could be a bridge to 77 collisions at a 

future e + e~ linear collider. ___ 

Over the past few years a comprehensive analysis of the heavy quark and 
I^IVI^IAH fill 1 11 *M 

vector meson 1 '"' 1 ' 1 ' production in UPC's has been made considering differ- 
ent theo reti cal approaches. For a recent review on heavy quark production in UPC's, 
see Ref. In particular, much e ffort h as been devoted to obtain signatures of the 
QCD Pomeron in such processes I ^ I Q I which can be used to constrain the QCD 
dynamics at high energies. On the other hand, recently the STAR Collaboration 
released the first data on the cross section of the coherent p production in gold - 
gold UPC's at y/s = 130 GeV El which provides the first opportunity to check 
the basic features and main approximations of the distinct approaches describing 
nuclear vector meson photoproduction. 

Here, we analyze the possibility of using UPC's as a photon-photon/nucleus 
collider, and study the production of vector mesons considering different QCD dy- 
namics. Probably, only the next generation of e + e~/ep/eA colliders will allow to 
discriminate between the distinct approaches for the QCD dynamics at high ener- 
gies. However, until these colliders become reality we need to consider alternative 
searches in the current and/or scheduled accelerators which allow us to constraint 
the QCD dynamics. The first analysis, summarized in Sec. [21 concerns to the pho- 
tonuclcar meson production, where we consider the QCD color dipole approach 
to describe its production in photon- nucleus process. That approach is suitable at 
high energies and allows take into account the corrections of partons saturation 
phenomenon (For recent reviews see, e.g. Ref. and nuclear effects in a simple 
and intuitive way. The basic quantity is the dipole cross section, which encodes all 
information on the interaction of the color dipoles with the nucleus target. The 
scattering of dipoles off nuclei is modeled through the Glauber-Gribov formalism. 
The conversion of dipoles on mesons is accounted by the meson wavefunction. Pre- 
dictions are given for the ion modes and energies of RHIC and LHC. A comparison 
with the recent STAR data on coherent p production is also presented. The second 
analysis, summarized in Sec. 01 treats the double meson production in two-photon 
processes, in particular double J/^>, p— J/ 1 ^ and double p production. In the case of 
heavy mesons, we consider the underlying dynamics at high energies as given by the 
perturbative QCD Pomeron, with the evolution described by the BFKL equation 
L^l Then we investigate the energy behavior in LO and NLO accuracy. The hard 
scale for the process is guarantee by the large mass of the heavy meson. For the 
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mixed light-heavy meson case, we use the QCD double logarithmic approximation, 
which depends on the gluon content of the light meson. The sensitivity to differ- 
ent gluon distributions is addressed. As a byproduct, we obtain the double p cross 
section using the previous results through the Pomeron-exchange factorization the- 
orem. In the last section we discuss the main features of the theoretical estimates 
and also address the backgroud processes. 



2. Photonuclear vector meson production at UPC's 

In heavy ion collisions the large number of photons coming from one of the colliding 
nuclei will allow to study photoproduction, with the photonuclear cross sections 
given by the convolution between the photon flux from one of the nuclei and the 
cross section for the scattering photon-nuclei. The final expression for the production 
of vector mesons in UPC's is then given by, 



& AA—>AAV 



(v 7 ^) = / du o iA ^va (W* A = 2ujy/S^j (1) 



where u> is the photon energy (u) m i n = m v /47j,m p ), ray is the meson mass and 
■sj Snn is the ion-ion c.m.s energy. For instance, the Lorentz factor for LHC is jl = 
2930, giving the maximum c.m.s. 'yN energy W 1 a ~ 950 GcV. In this process we 
have that the nuclei are not disrupted and the final state consists solely of the 
two nuclei and the vector meson decay products. Consequently, we have that the 
final state is characterized by a small number of centrally produced particles, with 
rapidity gaps separating the central final state from both beams. Moreover, due to 
the coherence requirement, the transverse momentum is limited to be smaller than 
Pt = \/2/Ra, where Ra is the nuclear radius. Therefore, these reactions can be 
studied experimentally by selecting events with low multiplicity and small total pr- 
The photon flux is given by the Weizsacker- Williams method The flux from 
a charge Z nucleus a distance b away is 

1 



duicPb ir 2 uib 2 



where is the Lorentz boost of a single beam and rj = ub/jL] Ko,i(v) are the 
modified Bessel functions. The requirement that photoproduction is not accompa- 
nied by hadronic interaction (ultraperipheral collision) can be done by restricting 
the impact parameter b to be larger than twice the nuclear radius, Ra — 1.2 A 1 ' 3 
fm. Therefore, the total photon flux interacting with the target nucleus is given by 
Eq. J2J integrated over the transverse area of the target for all impact parameters 
subject to the constraint that the two nuclei do not interact hadronically. An an- 
alytic approximation for AA collisions can be obtained using as integration limit 
b > 2 R Al producing 



dN (u) 2Z 2 a em [ rj 2 , 2 ^2/-^ 

(V) Kt [i]) + — [Ki (77) - K (77)) 



dcu it uj 



(3) 
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where fj = 2u> Ra/jl- 

Analyzing Eq. JQ) we obtain that the main input in the calculations of the 
vector mesons production cross sections in UPC's is the photonuclear cross section 
a (p/A — > V A) which is given by 

a ( 7 A - va) = P m A™(y = °)] 2 (1 + ^ r dt imi 2 (4) 

with train = (my/2w) 2 . Here, the scattering amplitude for meson production on 
photon- nuclei collisions at zero momentum transfer is labeled by A n uc- The quantity 
j3 is the ratio between the imaginary and real part of the amplitude. Using the 
analytical approximation of the Woods-Saxon distribution as a hard sphere, with 
radius Ra, convoluted with a Yukawa potential with range a = 0.7 fm, we obtain 
that the nuclear form factor reads as 1*^4, 

F{q = VW\) = -p? [sm(qR A ) - qR A cos{qR A )} 
Aq^ 

where po = 0.16 fm~ 3 . 

Currently, there are in literature different models for the photonuclear vector 
meson production cross section which differ basically in the treatment for the pho- 
ton and for its interaction with the target. The main aspect is that real photons have 
a complicated nature. In a first approximation, the photon is a point-like particle, 
although in field theory it may fluctuate also into a fermion pair (See discussions in 
Rcfs.LLLLJJ). i n the case where there is a photon transition in a colorless antiquark- 
quark pair, the propagation of this colorless hadronic wave packet in a nuclear 
medium can be treated either in the hadronic basis as a result of Gribov's inelastic 
corrections or in QCD in terms of the partonic basis, which are complementary. Lets 
briefly discuss these two representations (For a detailed discussion see Ref. 18 ). The 
time scale characterizing the evolution of a qq wave packet can be estimated based 
on the uncertainty principle and Lorentz dilation. The lifetime of the photon fluctu- 
ation is given by t c = vj (Q 2 + my), where v is the photon energy, my is the mass of 
the fluctuation and Q 2 is the photon virtuality. It is usually called coherence time. 
Using light-cone kinematics we can define the coherence length, which is given by 
l c = t c . Moreover, one cannot decide whether a ground state V is produced or the 
next-excited state V' , unless the process lasts longer than the inverse mass difference 
between these states. In the rest frame of the nucleus, this formation time is Lorentz 
dilated and is given by tf — 2v/(niy, — m v ). Similarly, we can define a formation 
length given by If =tf. In the hadronic basis, the same process looks quite differ- 
ent. The incident photon may produce different states on a bound nucleon, the V 
meson ground state or an excited state. Those states propagate through the nucleus 
experiencing multiple-diagonal and off-diagonal diffractive interactions, and even- 
tually the ground state is detected. According the quark-hadron duality, we expect 
that these two descriptions to be equivalent. However, as these two approaches have 
been used assuming different approximations, their comparison may provide a scale 



1 



1 + a 2 q 2 



(5) 
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for the theoretical uncertainty involved. Furthermore, it is important to emphasize 
that at high photon energy v, both l c and If greatly exceed the nuclear radius Ra, 
which implies in the partonic basis that the transverse size of the qq pair do not 
change during the interaction with the target. This enables one to introduce the 
QCD dipole picture ^11 where the process is factorized into the photon fluctuation 
in a qq pair and the dipole cross section. These aspects become the interpretation 
in the partonic basis more intuitive and straightforward than in the hadronic basis. 

Let us consider the scattering process jA — » V A in the QCD dipole approach, 
where V stands for both light and heavy mesons. The scattering process can be seen 
in the target rest frame as a succession in time of three factorizable subprocesses: 
i) the photon fluctuates in a quark-antiquark pair (the dipole), ii) this color dipole 
interacts with the target and, iii) the pair converts into vector meson final state. 
Using as kinematic variables the j*A c.m.s. energy squared s = W 2 , A = (p + q) 2 , 
where p and q are the target and the photon momenta, respectively, the photon 
virtuality squared Q 2 = —q 2 and the scaling variable x = (Q 2 + 4rnj)/(W^ A + Q 2 ), 
the corresponding imaginary part of the amplitude at zero momentum transfer reads 

as Eninr 

lmA( 1 A^VA)=J2 ! dzd 2 r^ h - h (z, r, QV^(£,r)*£* R (*, r) , (6) 
h,h 

where ^ r (z, r) and ^^(z, r) are the light-cone wavefunctions of the photon 
and vector meson, respectively. The quark and antiquark hclicities are labeled by h 
and h and reference to the meson and photon hclicities arc implicitly understood. 
The variable r defines the relative transverse separation of the pair (dipole) and z 
(1 — z) is the longitudinal momentum fractions of the quark (antiquark). The basic 
blocks are the photon wavefunction, , f 7 , the meson wavefunction, ^>)f L , and the 
dipole-target cross section, er^ sct - 

In the dipole formalism, the light-cone wavefunctions ^> h j l (z 1 r) in the mixed 
representation {z, r) are obtained through two dimensional Fourier transform of the 
momentum space light-cone wavefunctions, ^ h j l (z, k), which can be completely 
determined using light cone perturbation theory. On the other hand, for vector 
mesons, the light-cone wavefunctions are not known in a systematic way and they 
are thus obtained through models (For a recent detailed discussion see Ref. LLal). 
Here, we follows the analytically simple DGKP approach 1^21 w hich assumes that 
the dependencies on r and z of the wavefunction a re fact orised, with a Gaussian 
dependence on r (For a detailed discussion see 

Refs |22|24|) The 

main shortcoming 

of this approach is that it breaks the rotational invariance between transverse and 
longitudinally polarized vector mesons^! However, as it describes reasonably the 
HERA data for vector meson production, as pointed out in Ref. 1^1, we will use it in 
our phenomenological analysis. Finally, the imaginary part of the forward amplitude 
can be obtained by putting the expressions for photon and vector meson (DGKP) 
wavefunctions into Eq. ©. Moreover, summation over the quark/ antiquark hclic- 
ities and an average over the transverse polarization states of the photon should 
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be taken into account. The transverse component (the lon gitudinal one does not 
contribute for photoproduction) is then written as E0I3E322 



lmA T (s, t = 0)= d 2 r dza 1 ^ fy fr{z) exp 



+ (1 - z^K^er) + -±K {sr) ^(2,r) (7) 
my my r 



with cr^p get being the dipole-nucleus cross section. In the photoproduction case, 
£ = m y , where mf is the quark mass of flavour /. The corresponding parameters 
for the vector mesons wavefunctions (my, wr, fv, etc) are presented in Table 1 of 
Ref. El Following Ref. E2 

we have estimated contribution from real part for the 
photoproduction of vector mesons. This correction is about 3% for light mesons and 
it reaches 13% for J/ty at high energies Additionally for heavy mesons we have 
take into account the skewedness effects, associated to off-forward features of the 
process (different transverse momenta of the exchanged gluons in the t-channel), 
which are increasingly important in this case. 

In order to obtain the dipole-nucleus cross section we wil l assume the validity 
of the Glauber-Gribov picture ESI w hich allows to write EBE3 



^ Cleus (5- r 2 - A) 



d 2 b I 1 



exp 



^T A (b)^(x, r 2 ) 



(8) 



where b is the impact parameter of the center of the dipole relative to the center 
of the nucleus and the integrand gives the total dipole-nucleus cross section for 
a fixed impact parameter. The nuclear profile function is labeled by Ta(6), which 
will be obtained from a 3-parameter Fermi distribution for the nuclear density. The 
above equation sums up all the multiple elastic rescattering diagrams of the qq 
pair and is justified for large coherence length, where the transverse separation r of 
partons in the multiparton Fock state of the photon becomes as good a conserved 
quantity as the angular momentum, i. e. the size of the pair r becomes eigenvalue 
of the scattering matrix. Fo r th e dipole-proton cross section we follow the quite 
successful saturation model which interpolates between the small and large 
dipole configurations, providing color transparency behavior, <r^° on ~ 



, as r — > 

ctoi a t large dipole separations. The parameters 
of the model have been obtained from an adjustment to small x HERA data. The 
parameterization for the dipole cross section takes the eikonal-like forml^HI 



and constant behavior, c^ z r ° ton 



proton 
' dip 



(x, r ) = er 



1 — exp 



?sat(5) 



where the saturation scale Q 2 at defines the onset of the saturation phenomenon, 
which depends on energy. An additional parameter is the effective light quark mass, 
mf = 0.14 GeV, which plays the role of a regulator for the photoproduction (Q 2 = 0) 
cross section. The charm quark mass is considered to be m c = 1.5 GeV. 
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Table 1. The integrated cross section for vector mesons photonuclear produc- 
tion at UPC's at RHIC and LHC energies. 





HEAVY ION 


J/* (3097) 


0(1019) 


u (782) 


p(770) 


RHIC 


SiSi 


3.42 i^b 


612 fjb 


764 fib 


6.74 mb 




AuAu 


476 fjb 


79 mb 


100 mb 


876 mb 


LHC 


CaCa 


436 fjb 


12 mb 


14 mb 


128 mb 




PbPb 


41.5 mb 


998 mb 


1131 mb 


10069 mb 



Let's present the estimates from the QCD saturation model for vector meson 
photoproduction in the kinematical range of the colliders RHIC and LHC. It is 
noticeable that the energy dependence of the dipole-nucleus cross section [Eq. 
is strongly connected with the saturation scale Q1a(Wia) x ^ Qsnt- Namely, 
the saturation effects are larger whether the momentum scale is of order or larger 
than the correspondent size of the vector meson and the energy growth of the cross 
section is then slowed down. In Table ^ one shows the corresponding integrated 
cross section. For LHC energy, one considers lead (Pb) and calcium (Ca), whereas 
for RHIC one takes gold (Au) and silicon (Si). The integrated cross sections can 
be contrasted wit h t he theoretical estimations using GVDM plus Glauber-Gribov 
approach of Refs^^ as we n as the estimation of Ref. EH, which considers VDM 
plus a classical mechanical calculation for nuclear scattering and uses as input for 
the 7 p — > Vp reaction an extrapolation of the experimental DESY-HERA fits for 
mes on p hotoproduction. Initially lets consider the latter approach (see Tab. Ill in 
Ref. "ty . At RHIC energy and Si nucleus, our results are about 20 % lower for p and 
w, whereas gives a larger <p cross section and almost the same J/vE' cross section. 
However, the situation changes for gold nucleus, where the present results are about 
50 % larger than the estimates in Ref. At LHC energy and for Ca nucleus, our 
results gives higher cross section by a factor of order 10 %, whereas for lead nucleus 
the factor reaches a factor 2 for light mesons and almost similar for the J/^S meson. 
Basically, the values are quite similar for light nuclei. For heaviest nuclei, the results 
overestimate those ones in Ref. Gl wnen one considers light mesons and become 
similar for the J/\l/ case. These features can be understood through the theoretical 
procedure when considering the nuclear scattering. 

iti • rm 

Now we compare our results ( with those ones in Ref. ^-H where the main focus 
is on the p and J /iff production. For p the predictions are computed only for RHIC 
energy ^Jsmn = 130 GeV and we will consider it later on. We can anticipate that 
their results are closer to ours since a Glauber-Gribov approach is used in describ- 
ing the scattering on nuclei. For J/vf the theoretical approach for the photonuclear 
production was the collinear QCD double logarithmic approximation, where the 
7A — > J/*f?A cross section is directly proportional to the squared nuclear gluon 
density distribution 1^4 There, it was considered an impulse approximation (no nu- 
clear shadowing) and a leading twist shadowing version. The impulse approximation 
gives a larger cross section at central rapidity (about a factor 4 higher for Ca and 
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A STAR DATA 

CUTS: K0.02 
lyl<3 



CUTS: K0.02 
lyKl 



100 150 200 
W NN [GeV] 



1 00 1 50 200 
W„ [GeV] 



Fig. 1. Energy dependence of coherent p meson production in gold-gold in UPC's at RHIC 
(y/ s NN = 130 GeV). Experimental data from STAR Collaboration. 



a factor 6 for Pb), while at fragmentation region both approximations match each 
other at LHC energy for Ca and Pb nuclei. Our results are closer to their impulse ap- 
proximation, which suggests nuclear shadowing could be weak for J/^> production. 
This feature can be easily tested in the first experimental measurements of coherent 
J/^f production on UPC's at LHC. Concerning the integrated cross section, they 
found 0.6 mb for Ca nucleus and 70 mb for Pb at LHC. Our results are 0.44 mb 
and 41.5 mb, respectively. Thus, our results are about 15 % lower for Ca and also 
40 % lower for Pb. The difference between the predictions comes from mostly from 
the distinct QCD approaches considered used and the different photon flux in the 
UPC calculation. 

Recently, the STAR Collaboration published the first measuremente for coherent 
p production in gold - gold UPC's at y/s = 130 GeVl^U providing the first opportu- 
nity to test the distinct approaches describing nuclear vector meson photoproduc- 
tion. Our predictions are presented in Fig. ^ As the parameters in our approach 
were constr ained i n the study of the photoproduction of vector mesons in jp(A) 
interactions l ^3 | 24| our predictions are parameter free. In the range of rapidities 
\y\ < 3, at energy ^snn = 130 GeV, we have found cr sa t(— 3 < y < 3) = 410 mb in 
good agreement with the STAR measurement ctstar(~3 < y < 3) = 370 ± 170 ±80 
mb. For the cut \y\ < 1, we have obtained cr sat (— 1 < y < 1) = 221 mb, whereas the 
STAR result is ct S tar(-1 < V < 1) = 140 ± 60 ± 30 mb. The values presented here 
are somewhat similar to the ones obtained in Ref Hil, which uses the generalized 
vector dominance model (GVDM) and the Glauber-Gribov approach, including in 
addition the finite coherence length effects. 
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3. Vector meson production in two-photon processes at UPC's 

Rclativistic heavy-ion collisions are a potentially prolific source of 77 collisions at 
high energy colliders. The advantage of using heavy ions is that the cross sections 
varies as Z 4 a 4 rather just as a 4 . Moreover, the maximum 77 collision energy W 11 
is 2j/Ra, about 6 GeV at RHIC and 200 GeV at LHC, where Ra is the nuclear 
radius and 7 is the center-of-mass system Lorentz factor of each ion. In particular, 
the LHC will have a significant energy and luminosity reach beyond LEP2, and 
could be a bridge to 77 collisions at a future e + e~ linear collider. For two-photon 
collisions, the cross section for the reaction AA — » AA Vi V2 will be given by 

— -m(uji) / — - n 2 (w 2 ) (J 11 -,v 1 v 2 {W = V4^iw 2 ) , (10) 

where the photon energy distribution n(u>) is calculated within the Weizsacker- 
Williams approximation In general, the total cross section AA — > AA77 — > 
AAX, where X is the system produced within the rapidity gap, factorizes into the 
photon-photon luminosity and the cross section of the 77 interaction, 



dr 

r2 



&aa^aav 1 v 2 ( s ) = / dr — — &~~_ tVl v a (s), (11) 



where r = s/s, s = W is the square of the center of mass (c.m.s.) system energy 
of the two photons and s of the ion-ion system. The 77 luminosity is given by the 
convolution of the photon fluxes from two ultrarelativistic nuclei: 

c££ 77 (r) f dx 



dr 



-f(x)f(r/x), (12) 



where the photon distribution function f(x) is related to the equivalent photon 
number n(uj) via f(x) = (E/oS)n(xE), with x = lu/E and E is the total energy 
of the initial particle in a given reference frame. The remaining quantity to be de- 
termined in order to proceed is the quantity f(x), which has been investigated by 
several groups (for more details, see e.g. ™). Here, we consider the photon distri- 
bution of Rcf. ~ providing a photon distribution which is not factorizable. The 
authors of tiJJ produced practical parametric expressions for the differential lumi- 
nosity by adjusting the theoretical results. The comparison with the complete form 
is consistent within a few percents. The approach given above excludes possible fi- 
nal state interactions of the produced particles with the colliding particles, allowing 
reliable calculations of UPC's. Therefore, in order to estimate the double vector 
meson production cross section in UPC's it is only necessary to consider a suit able 
QCD model for this process in 77 collisions (For a recent analysis see Ref . ^|) . 

Let's start our dicussion considering the double J/^S production. In Ref. 
the double J/ 1 ^ production in 77 collisions has been proposed as a probe of the 
hard QCD Pomeron. At the photon level we have that the calculation of the cross 
section can be made perturbatively, due to the presence of the charm mass. The 
non-perturbative content is provided only by the J /\& light-cone wave function, 
which is well constrained through the experimental measurement of its lcptonic 
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width T j/^^i+i- . In the following we use the high energy factorization and the 
BFKL dynamics in order to perform estimates for the referred reaction. Assuming 
a small t approximation, the total cross section for this process can be written as, 

1 da (7 7^ J/* J/*) 



Otct (7 7 -> J/® J/®) = 
with 



Bj/y J/9 dt 



(13) 



t=o 



It ~~ 16tT P \ Bj / xVJ l* ' ' ^ 14 ) 

where Bj/*,//* is the corresponding slope parameter, which we assume as being 
Bj/qi j/<s — Tn~ 2 (See discussion in Ref.1% In the general case, the imaginary part 
of the scattering amplitude is given as follows^!, 

ImA{W 1T t)=J {k + g/2)2 {k _ q/2)2 , (15) 

where W is the center of mass energy of the two photon system and the photon- 
meson impact factor is denoted by <& 7 ,7/^. At the Born level $ = $ and the reaction 
is described by the two-gluon exchange, which have transverse momenta q/2 ± k 
and where the momentum transfer is t = —q 2 . When considering the complete 
gluon ladder contribution, the quantity $ contains the impact factor and the gluon 
emission on the ladder, which is driven by the QCD dynamics. At the LLA level, 
the BFKL ladder contribution for the i-channel exchange provides the following 
expression for it, 

r k 2 

$ 7j/ *(W^,k,q) = J d 2 k' -^1C{W^,k,k' ,q)%j,*{k' ,q) , (16) 



where 



$ 7 j/*(fe,q) = C 



q 2 m 2 + k 2 



(17) 



and K(W* 7 ,k, fe', q) is solution of the LLA BFKL equation at t ^ 0. The Eq. (fTYjl 
defines the impact factor in the nonrelativistic approximation, where it is assumed 
that the quark and antiquark have the same longitudinal momentum. A comparison 
between the HERA data and the theoretical predictions for J/^> photoproduction 
using this approximation has found good agreement We have considered the 
following parameters for the further calculations: m c = = 1.55 GeV, C — 

y/a em a s {fJ, 2 ) e c 1 7r m c fj/y, with e c = 2/3 and fj/q, = 0.38 GeV. 

In the Born two-gluon level we have that JC = 5 2 (k — k) and the cross section 
is energy independent. On the other hand, in order to perform a LLA BFKL cal- 
culation, the following solution for the evolution equation in the forward case was 
considered El, 

1 1 (w 2 Y F '" 2 (fc 2 /.fc' 2 > 



^2^ak 2 k rz l\n{W 2 /s 
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where ujp = -2-2^-4 In 2 and a = 28 £(3). The Pomeron intercept is given by 
1 + ljp in the leading logarithmic approximation, which depends on a s . We have 
taken s = 1 GeV 2 . The results are sensitive to the choice for the intercept, providing 
an enhancement of the total cross section by one/two orders of magnitude in relation 
to the Born level at the considered energy range. The behavior presented by the 
approximation above is considerably steep on energy. Consequently, it is important 
to evaluate the corrections which comes from the next-to-leading order corrections 
to the BFKL equation or those associated to unitarity effects. In order to simulate 
the NLO effects, we have used a value top = 0.37. The effective behavior obtained 
is given by W^ , with A = 0.29, which value is somewhat close to A ~ 0.23 — 0.29 
obtained inl^. A more detailed analysis of the double J/^ production at the photon 
level can be found in Ref. El 

Having determined the 77 cross section we can estimate the double J/\I/ cross 
section in UPC's. In Ref. El the energy dependence of the cross section was investi- 
gated, considering the Pb + Pb collisions and different approximations for the 77 
subprocess. Shortly, the LLA BFKL solution predicts large values of cross section 
for LHC energies in comparison with the Born approximation El The ratio between 
the BFKL prediction and the two-gluon cross section is around 40 at yfs = 5500 
GeV. Taking into account non-leading corrections, this ratio is reduced to approx- 
imately 1.7. Therefore, these results indicate that if the NLO corrections to the 
BFKL approach are account in the approximation considered here, a future exper- 
imental analyzes of double J /i$> production in ultraperipheral heavy ion collisions 
could not easily constrain the QCD dynamics, since the BFKL(NLO) result is sim- 
ilar to the Born prediction. Only accurate measurements would allow discriminate 
between the two cases. Moreover, it is important to emphasize that the NLO correc- 
tions can be larger, implying a full suppression of energy enhancement associated 
with the iteration of gluons in the t-channel present in the QCD dynamics at high 
energies. However, if the energy dependence of the 77 cross section was driven for 
a large intercept, closer to the LO prediction, the analysis of UPC's can be useful. 
It is worth to mention that our results are quite similar to those ones computed for 
the FELIX proposal at LHC For PbPb collisions with energies of center of mass 
equal to y/s = 5500 A GeV, luminosities of Caa — 4.2 x I0 26 cto~ 2 s _1 are consid- 
ered. Consequently, during a standard 10 6 s/ month heavy ion run at the LHC, we 
predict approximately 30, 50 and 1100 events for Born, BFKL (MOD) and BFKL 
(LO), respectively El. It is interesting to compare these predictions with the results 
for double J "/\f r production at LEP2, since the 77 center of mass energies are simi- 
lar. Considering an integrated luminosity of 500 pb -1 in three years and \fs = 175 
GeV, almost 70 events are expected taking into account the non-leading corrections 
to the BFKL approach 1^21. Therefore, we predict a large number of events in ultra- 
peripheral heavy ion collisions, allowing future experimental analyses, even if the 
acceptance for the J/\l/ detection being low. It should be stressed that the LHC 
probably will operate in its heavy ion mode only four weeks per year. 
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Let's now discuss the pJ/^ production in UPC's . At the photon level, this 
process can be calculated similarly as made for the elastic J/\V photoproduction off 
the proton Basically, the differential cross section is given by^l 

^n^m = c ^ 3 9> : (f)? r e r [xGW)]2 e - W i (19) 

dt 3M ,//* 

where C denotes a product of corrections factors. The two-photon cms energy is 
denoted by W 77 , where x — Mj, 9 /W^ and Mj/m is the heavy meson mass. The 
mass scale is given by p? — M 2 ^/A. In the small-t approximation, the slope is 

estimated to be B pJ /^, = 5.5 ± 1.0 GeV" 2 E3. The light meson-photon coupling 
is denoted by g 2 = 0.454 and the heavy meson decay width into a lepton pair 
is r ee . Moreover, xG p (x,p 2 ) is the gluon content of the light meson, which can be 
contraint from experimental data for the cross section. Currently, there are different 
parameterizations for this distribution in the literature, associated to distinct QCD 
dynamics. In Ref. ^we investigate if an experimental analysis of the pJ/^ allows 
discriminate among distinct parameterizations for the gluon content of the light 
meson. In particular, we use both LO and NLO GRS parameterization Ell the 
SaSID parameterization ESI and a phenomenological Regge motivated ansatz EH. 
The latter is given by x (x, p 2 ) = x Q G p (x , p 2 ) (^)^ JF , with x = 0.1 and 

the Pomeron intercept is considered as ojp = 1.25, in agreement with the effective 
power in the HERA data. For the nuclear case we allow for a higher intercept, in 
order to simulate a BFKL-like behavior, motivated by the studies in double J/\& 
production 3 . The normalization, x G p ([xq, M 2 ^/4), is given by the NLO-GRS 
parameterization. The differences among the parameterizations are sizeable already 
at photon level as a consequence of distinct effective exponent A. The GRS LO 
one presents the steeper behavior, followed by SaSID (A = 0.3038). The Regge 
motivated and GRS NLO parameterizations are more close since the intercept for 
the Regge ansatz is similar to the effective power of GRS NLO (A ~ 0.228). The 
results for the pJ/^ production cross sections at the photon level can be found in 
Ref.El. 

For p J/^> meson production in UPC's at LHC energy range we have that an 
upper bound is obtained using the SaSID parameterization and the lower one by 
the Regge ansatz with the lower Pomeron intercept ujp — 1.25. At the LHC energies 
y/s = 5500 GeV, the cross section takes values between 680 nb and 2.6 mb, show- 
ing the process can be used to discriminate among parameterizations. It should be 
noted that the ultraperipheral cross section is dominated by not so high two-photons 
energies. Although the deviations among the models are quite visible at very high 
energies, in the nuclear cross section such deviations are less sizeable. This is directly 
related to the effective two-photon luminosity, which peaks at smaller W 77 . There- 
fore, correct estimates should include careful treatment of the low energy threshold 
effects. The nuclear cross section is enhanced in relation to the e + e~ case, which 
provides values of hundreds of pb's in contrast with units of mb in the peripheral 
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Table 2. The double p production cross sections in UPC's at 
LHC (y/s = 5500 GeV) for PbPb. 



SCENARIO 


o(AA -> AApp) (nb) 


BFKL(LO) + GRS(LO) 


25 x 10 3 


BFKL(LO) + GRS(NLO) 


4 x 10 3 


BFKL(LO) + Rcggo {u) P = 1.25) 


810 


BFKL(MOD) + GRS(LO) 


125 x 10 4 


BFKL(MOD) + GRS(NLO) 


20 x 10 4 


BFKL(MOD) + Regge (ui F = 1.25) 


405 x 10 2 



nuclear case. Comparing our estimates with the FELIX parameterization it was 
found the values predicted by FELIX collaboration are ever larger than obtained 
using the approach presented here, independently of the value of slope used. As a 
cross check, we also calculated the cross section for pJ/^ production in UPC's using 
the same kinematical cut proposed in Ref. ESI In such a situation, the agreement 
was much better. 

As a byproduct, using the Pomeron-exchange factorization theorem we can 
use the results for double J/ 1 ^ production and for p J/^ production to obtain a rough 
estimate of the double p production in ultraperipheral heavy ion collisions. In prin- 
ciple, this theorem cannot be applied for BFKL Pomeron, since it is not an isolated 
Regge pole. As pointed out in Ref.l^ (See also Ref. ^j), only when the asymptotic 
freedom is incorporated into the BFKL equation is that the BFKL Pomeron can be 
expressed in terms of a series of isolated poles in the angular-momentum plane, with 
the contribution of each isolated pole satisfying the factorization theorem. An open 
question is how much the factorization is violated in different processes. In what 
follows we estimate the double p cross section assuming the Pomeron-exchange fac- 
torization theorem. These predictions should be compared with those obtained in 
a QCD calculation. Preliminary results presented in Ref. ESI shown that the viola- 
tion is not large. However, more detailed studies are necessary in order to quantify 
this violation. Given the assumption that single Pomeron exchange dominates, the 
following relation among the total cross sections can be posed 

■KT,~»)- 'y * yaffil - m 

Using the factorization theorem for the slopes, we estimate B pp w 12.1 GeV -2 . In 
the Table[21we present our predictions for the double p production in ultraperipheral 
AA collisions, considering different scenarios for the QCD dynamics. We can see 
that there is a large range of possible values for the cross section, which implies 
that future experimental data are essential to constraint the dynamics as well as 
the violations of the factorization theorem. It is important to emphasize that our 
predictions agree with those presented in Ref. ESI 
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4. Summary 

Before presenting a summary of the main results, lets discuss the background pro- 
cesses and experimental separation. In Ref.^a detailed analysis of the experimental 
separation between photoproduction and two-photon interactions was presented. 
There, the authors have estimated that the two-photon cross sections are at least 
1000 smaller than the photoproduction cross section, which makes the experimen- 
tal separation between the two interactions very hard. Our calculations indicate 
that the inclusion of the QCD Pomeron effects implies higher cross sections at two- 
photon level and, consequently, larger cross sections in ultraperipheral collisions. 
Therefore, the inclusion of these effects implies that, in general, the contribution of 
two-photon interactions is non-negligible. However, the experimental separation of 
two-photon process still remains a challenge. In principle, an analysis of the impact 
parameter dependence should allow to separate between the two classes of reactions, 
since two-photon interactions can occur at a significant distance from both nuclei, 
while a photonuclear interaction must occur inside or very near a nucleus. 

Other possible background process is that associated to Pomeron-Pomeron in- 
teractions. It has been verified that such reactions would be non-negligible for light 
ions, while they are significantly suppressed for heavy ions 1^21. In the particular 
case of the double heavy meson production this contribution deserves more detailed 
studies, since the current treatments rely on the Regge formalism instead of a QCD 
approach. An additional contribution in two photon ultraperipheral collisions is the 
meson production accompanied by mutual Coulomb dissociation. It has been esti- 
mated that these reactions increase the total cross section for an amount of ~ 10% 
at RHIC/LHC energies ^1 We disregarded this contribution in the present calcu- 
lations, since that it would be smaller than the associated theoretical uncertainties. 

Let's consider the experimental feasibility of photonuclear production. Although 
the vector meson photoproduction at AA collisions to be a small fraction in compari- 
son to the total hadronic cross sections, the experimental separation of these reaction 
channels is possible. As photoproduction is an exclusive reaction, A+A — > A+A+F, 
the separation of the signal from hadronic background would be relatively clear. 
Namely, the characteristic features in photoproduction at UPC's are low px meson 
spectra and a double rapidity gap pattern. Moreover, the detection (Roman pots) 
of the scattered nuclei can be an additional useful feature. In hadroproduction, the 
spectra on transverse momentum of mesons are often peaked around meson mass, 
Pt ~ my. An experimental cut px < 1 GeV would eliminate most part of the 
hadronic background. Hence, the rapidity cut would enter as an auxiliary separa- 
tion mechanism. This procedure is specially powerful, since there will be rapidity 
gaps on both sides of the produced meson. 

As a summary, the cross sections for the A + A — > A + A + V (V = p,u>,<j>, J/^f) 
process were computed and theoretical estimates for scattering on both light and 
heavy nuclei are given for RHIC and LHC energies. The rates are very high, mostly 
for light mesons and at LHC energies. As an important result, we compare our 



Vector Meson Production in Ultraperipheral Heavy Ion Collisions 15 



prediction for the coherent p meson production with RHIC data at ^snn = 130 
GeV. The corresponding results are in good agreement with the experimental re- 
sults when considered the cuts on momentum transfer and on rapidity. Moreover, 
we have computed the double vector production in UPC where one considers the 
perturbative QCD Pomeron, described by the BFKL equation. For the mixed light- 
heavy meson production, we use the DDLA pQCD approximation, which depends 
on the gluon content of the light meson and obtain the double p cross section using 
the previous results through the Pomeron-exchange factorization theorem. 

As a final remark, in a similar way the strong electromagnetic fields generated 
by high-energy protons allow us to study photon - nucleon processes in proton - 
proton interactions in a large kinematical range. These events can be experimen- 
tally studied by selecting events with low multiplicity and small total transverse 
momentum. We have studied this case for the photoproduction of vetor meson s in 
proton- (anti)proton in Ref. ^and the photoproduction of heavy quarks in Ref. 
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